Disordered contact process with asymmetric spreading 
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An asymmetric variant of the contact process where the activity spreads with different and in- 
dependent random rates to the left and to the right is introduced. A real space renormalization 
scheme is formulated for model by means of which it is shown that the local asymmetry of spread- 
ing is irrelevant on large scales if the model is globally (statistically) symmetric. Otherwise, in the 
presence of a global bias in either direction, the renormalization method predicts two distinct phase 
transitions, which are related to the spreading of activity in and against the direction of the bias. 
The latter is found to be described by an infinite randomness fixed point while the former is not. 
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I. INTRODUCTION 

The contact process [l|, Q is a simple of model of 
an epidemic that is defined on a lattice of binary vari- 
ables which can be either inactive (healthy) or active 
(infected). The dynamics consist of two kinds of com- 
peting moves: infection of adjacent healthy sites by in- 
fected ones and spontaneous healing of the latter. Be- 
sides being a starting point for designing more realis- 
tic models of epidemics it is a thoroughly studied (but 
non-soluble), paradigmatic model of systems undergo- 
ing a non-equilibrium phase transition from a fluctuating 
phase to an absorbing one 0-01 • In. the case of transla- 
tional invariance, the phase transition belongs to the uni- 
versality class of directed percolation 3, although this 
type of critical behavior is rarely observed in real systems 
[6| since, according to the Harris criterion, it is unstable 
against weak quenched disorder in dimensions d < 4 0. 
The understanding of the behavior of the model in the 
presence of quenched disorder is therefore of great impor- 
tance. Early Monte Carlo simulations showed, in agree- 
ment with phcnomenological considerations, that the dis- 
order leads to anomalous dynamics outside of the critical 
point 041] and, at the critical point, the different ob- 
servables scale with some power of the logarithm of time 
rather than the time itself Q. Below the critical point, 
namely, the density of active sites and the survival proba- 
bility show a power-law decay with dynamical exponents 
varying continuously with the control parameter, while, 
above the critical point, the size of the set of active sites 
is growing sub-linearly in one dimension when the pro- 
cess was started from a single active seed [l3|- These 
phases are analogues of the disordered and ordered Grif- 
fiths phases of magnetic systems respectively, and 
the slow dynamics are related to the occurrence of so 



called rare regions which are locally in the opposite phase 
with respect to the majority of the system. In the sub- 
critical Griffiths phase, the rare regions are locally super- 
critical clusters and, although, being exponentially rare 
(in their size) , due to their exponentially large extinction 
time they are able to change the usual exponential tem- 
poral decay of the density to an algebraic one. Instead, 
in the super-critical Griffiths phase, the creeping motion 
of the front of the active cluster is caused by rare sub- 
critical regions which impede the spreading of activity 
for long times. A substantial progress in the quantitative 
description of the critical point was the adaptation of a 
strong disorder renormalization group (SDRG) scheme 



12 1 to the model by Hooyberghs, Igloi and Vanderzande 
13| . The formal description of the model by the SDRG is 
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essentially identical with that of the random transverse- 
field Ising model [l^ and, in one dimension, yields loga- 
rithmic critical dynamics and the complete set of critical 
exponents which are universal, i.e. independent of the 
form of disorder. These predictions were confirmed by 
Monte Carlo simulations and found to be valid also for 
relatively weak disorder, where the initial steps of the 
SDRG procedure are highly approximative [T5| . 

In the SDRG treatment of the model, the activity is 
assumed to spread through a given link in both direc- 
tions with the same rate (that varies from link to link), 
while in the numerical simulations, slightly differently, 
the infection from a given site occurs in both directions 
with equal rates. So, in the latter case, the spreading 
of activity through a link is non-symmetric (or biased) 
and one may ask whether this case can be treated by 
an appropriately generalized SDRG scheme. More gen- 
erally, one could pose the same question as well as the 
question what is the nature of the phase transition in 
an asymmetric variant of the disordered contact process 
where the infections rates \ij assigned to a directed link 
(i,j) are completely uncorrelated random variables. In 
case of a translationally invariant model with a bias (i.e. 
different infection rates in the two directions) the phase 
transition is known to belong to the directed percolation 
universality class as in the unbiased model |16| . But, the 
breaking of the translation symmetry, for instance, by 
putting an active wall to the system, leads to differences 
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compared to the unbiased case, see Ref. [17|. So, one ex- 
pects that that in the simultaneous presence of a global 
bias and quenched disorder the nature of phase transition 
is different from that of the homogeneous, biased contact 
process. 

The aim of the present work is to study the disor- 
dered, asymmetric contact process in one dimension by 
working out and applying an SDRG scheme. We will 
show that, if the model is statistically symmetric, then 
the local asymmetry of links is irrelevant and the critical 
behavior is identical to that of the symmetric model. If, 
however, the model is globally asymmetric, i.e. there is 
a bias in either direction, two distinct phase transitions 
arise, which are related to the spreading of activity in and 
against the direction of the bias. These results will be ex- 
pounded in the rest of the paper, which is organized as 
follows. In Section ini the precise definition of the model 
is given and the SDRG scheme is introduced. The case 
of weak asymmetry is analyzed in Section IIII[ while the 
behavior of the globally symmetric and the biased model 
is discussed in Section |IV] and Section |Vl respectively. 
The results are summarized in Section [VII and some cal- 
culations are presented in the Appendix. 
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FIG. 1: Left: Illustration of the merging of two sites con- 
nected by a link with large infection rates. Right: Schematic 
spectrum of the rate matrix of the two-site block. 

the symmetric model, there are two types of reduction 
steps, and this structure of the SDRG with appropriate 
modifications will be kept also for the asymmetric contact 
process. 

A. Cluster merging 



II. THE ASYMMETRIC CONTACT PROCESS 
AND THE SDRG SCHEME 

Let us consider a one-dimensional lattice where each 
site can be in two states: either active (infected) or 
inactive (healthy). The asymmetric contact process is 
a continuous-time Markov process on this state space, 
which consists of the following transitions occurring in- 
dependently. Active sites become inactive with a rate /i; 
on site i, and, if site i is active, it infects the neighbor- 
ing site on its right with a rate Xi and the neighboring 
site on its left with a rate Ki-i. The rates and Ki 
will be termed as forward and backward infection rates, 
respectively. The transition rates {fj,i, Xi, Hi} are time- 
independent random variables, the distribution of which 
will be specified later. 

The SDRG transformation of the disordered, symmet- 
ric contact process is a real space renormalization 
transformation which sequentially eliminates the quickly 
relaxing degrees of freedom and replaces the original sys- 
tem by a reduced one that, however, preserves the slowly 
relaxing degrees of freedom of the original system. For- 
mally, finite blocks of sites are considered, the time evo- 
lution of which are governed by a master equation 



dtP = PQ, 



(1) 



where the probabilities of states are arranged in the row 
vector P{t) = {pi{t),p2{t), . . . ) and Q is the rate matrix 
(or infinitesimal generator) of the process. In the SDRG 
procedure, the blocks are selected in which the spectrum 
of Q shows a separation of niveous. Then the higher- 
lying states are discarded and the block is replaced by 
a simpler one having the same low-lying niveous. For 



First, let us consider an isolated block of two sites 
shown in Fig. [TJ and assume that the infection rates in 
both directions are much larger than the recovery rates, 
Ai,Ki ^ /ii,/i2- In this case the two sites are most of 
the time either both active or both inactive. This sug- 
gests that the block can be substituted by a single giant 
site (or cluster) with an appropriate recovery rate fl. In- 
deed, analyzing the spectrum of the rate matrix of the 
block, see the Appendix, one obtains that the two higher- 
lying niveous are well separated from the two lower-lying 
ones. Note that at least one of the eigenvalues is al- 
ways zero since the rate matrix is a stochastic matrix. 
For an isolated single site with recovery rate fl, the only 
non-trivial eigenvalue is —fl. So, the effective recovery 
rate of the giant site is identified with the magnitude of 
lowest non- trivial eigenvalue of the rate matrix or, equiv- 
alently the (in magnitude) smallest root of Eq. IA2I If 
Ai, Ki ^ A'i!M2i one obtains the approximate form: 



M - 



where the combined infection rate tJi is defined as 

AjKj 

= — . 

Xi + Ki 



(2) 



(3) 



Note that if the infection rates Xi and Ki differ by or- 
ders of magnitudes then Ui gives roughly the smaller one 
of them. We can see that this type of reduction of the 
two-site block to a single site is justified only if both in- 
fection rates on the link (or, equivalently uJi) are much 
larger than the recovery rates [23, and the life-time of 
activity in the two-site block is controlled by the small- 
est infection rate. Consequently, in a totally asymmetric 
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FIG. 2: Illustration of the elimination of a site with a large 
recovery rate /^2- 
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FIG. 3: Schematic spectrum of the rate matrix of a three- 
site block in the case of complete separation (left) and partial 
separation (right). 



contact process where infection spreads unidirectionally, 
giant clusters with low decay rate never form, hence no 
rare region effects are expected to emerge (unless the ini- 
tial recovery rates are not bounded away from zero). 

Besides the transition rates, we also keep track of the 
mass of clusters, i.e. the number n of original sites com- 
prised by a giant site. Obviously, this quantity trans- 
forms as 

ri = 71l + 712. (4) 

B. Cluster elimination 

Next, let us consider an isolated three-site block shown 
in Fig. [21 where the recovery rates on the two lateral sites 
are zero. If the recovery rate [ii is much larger than the 
infection rates. 



and obtain the approximate expressions 



[ll > Ai,Ki,A2,K2, 



(5) 



the middle site will be most of the time inactive and 
transmits infection only during its rare flashes of activity. 
It is thus reasonable to eliminate this site and connect the 
lateral sites by a direct link with some effective (reduced) 
infection rates. Solving the eigenvalue problem of the 
rate matrix of the block, see the Appendix and Fig. [3l 
one sees that the higher-lying 4 niveous, being in the 
order of ^2, are well separated from the lower- lying 4 
niveous among which two niveous have zero eigenvalue. 
This lower-lying bunch of niveous is to be compared with 
the spectrum of a two-site block with infection rates A, k 
and zero recovery rates. The eigenvalues of the latter 
are (0, 0,— A, — k), so the effective infection rates of the 
three-site block are identified as the magnitude of the 
two lowest-lying non-trivial eigenvalues. If the condition 
in Eq. ([5]) is fulfilled, we have 



A = SA 1- Wl 



Ai A2 



/i2 + Ai + A2 

SA = ^ (6) 



K1K2 



+ Ki -t- K2 



(7) 



A1A2 



(8) 



while the other non-trivial eigenvalues are 0(^2) ■ Here, 
we speak of a complete separation of niveous. As we 
will see later, the cluster elimination step will be ap- 
plied in the SDRG procedure only if ^2 > wi,a;2 and, 
due to the broad distribution of rates we have typically 
LUi ~ min{Ai, k^}. If the two links of the block are biased 
in opposite directions, e.g. Ai, K2 ^ /i2 and ki, A2 > /i2, 
then the spectrum still displays a weaker form of com- 
plete separation (30l |. But if the two links are biased in 
the same direction, e.g. to the right. 



Al,A2>/^2, 'tl, <C /i2, 



(9) 



then only the lowest (non-trivial) eigenvalue separates 
from the higher ones. The corresponding rate k ~ 
KiK2/lJ'2 can still be interpreted as an effective spread- 
ing rate of the activity from right to left but by this 
reduction, the information on how the activity spreads 
from left to right is lost. We call this situation a partial 
separation of levels. 

We will also keep track of evolution of the number I 
of original links which are incorporated into the effective 
links. This quantity is dual to the mass of giant sites and 
transforms in the cluster elimination step as 



l=ll+l2 



The SDRG scheme 



(10) 



After introducing the elementary reduction steps, we 
are in a position to formulate a SDRG scheme for the 
asymmetric contact process. First, the largest one among 
the rates {wi, fii}, which will be denoted by il, is selected. 
If is a combined infection rate (recovery rate) then 
the cluster merging step (cluster elimination step) is per- 
formed. The above steps are then iteratively applied to 
the rcnormalized system. This recursive procedure grad- 
ually reduces the cut-off rate fl, which amounts to that. 
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roughly speaking, the fast processes occurring on a time 
scale shorter than are eliminated. The basic dynam- 
ical relation between the time and the length scale can be 
inferred from the asymptotic dependence of the cut-off 51 
on the fraction of active (non-decimated) degrees of free- 
dom [l2|. Note that, if the model is symmetric (Ai = k;) 
the above SDRG procedure reduces to the that applied 
inRcf. 

Before we proceed with the analysis of the SDRG 
transformation, we need to recapitulate some of the re- 
sults of the symmetric contact process obtained by the 
SDRG technique. For further details we refer the reader 
to Ref . . In the case of the symmetric model, the dis- 
tribution of the logarithm of rates is broadening without 
limits under the SDRG transformation, therefore the con- 
ditions of the approximative reduction steps are fulfilled 
better and better and the procedure becomes asymptot- 
ically exact. Due to that the rates at different sites re- 
main independent provided they were so initially, it is 
sufficient to deal with the evolution of the probability of 
rates under the renormalization and the problem is an- 
alytically tractable. The asymptotic transformation of 
rates in Eqs. ^ and ([5]) is identical to that of the cou- 
plings and external fields of the random transverse-field 
Ising chain up to a constant factor which is irrelevant at 
the critical point [l^. The asymptotic solution of the 
renormalization equations of the latter problem at the 
critical point, which is at the self-dual point of the trans- 
formation, have been found in Ref. |14| . Here, in the so 
called infinite randomness fixed-point of the transforma- 
tion, both the logarithm of infection and recovery rates 
scale as 



of infection rates relative to that of recovery rates by 



|lnA|^ lln^l^e 



1/2 



(11) 



and the mass of clusters as well as the length of effective 
links grow as 



(12) 



where df = {I + \/5)/4 — 0.809 • • • is the fractal dimen- 
sion of clusters. The Griffiths phase which is located on 
the sub-critical side of the critical point is characterized 
by a line of fixed points where the rates scale as 



llnAl^e, A^r% 
while the mass and the length behave as 
n In^, I ^ ^. 



(13) 



(14) 



The dynamical exponent z appearing here depends on 
the distribution of rates of the original model [l^. In 
the super-critical Griffiths phase A and n, as well as n 
and I are interchanged in the above relations. 



III. WEAK ASYMMETRY 

Returning to the asymmetric contact process, we intro- 
duce a control parameter which measures the magnitude 



A = In uji — In /ii 



(15) 



where the over-bar denotes an average over sites, and a 
local asymmetry parameter as 



Qi — In Ai — In Ki 



(16) 



First, let us consider a weak, asymmetric perturbation of 
the symmetric model: the parameters arc uncorrclatcd 
random variables written in the form 



Icti, 



(17) 



where / is an infinitesimally small global constant and 
ai is 0(1). Let us assume, furthermore, that the trans- 
formation is close to the fixed-point, so that the ap- 
proximate but asymptotically exact rules in Eqs. ([2]) 
and ([8]) can be used. Now, it is expedient to use the 
variables (idi^ai) instead of (Ai,Ki) in terms of which 
the old variables are expressed as Aj = LJi{l + e"') and 
Ki = uJi{l + e~°*). In a cluster elimination step these 
variables transform (asymptotically) as 



a ~ ai + a2 
w ~ /(ai,a2) 



UJ1UJ2 



f{ai,a2) = 1 + 



1 + e''i+'^2 ' 



(18) 
(19) 

(20) 



while in a cluster merging step they still follow the rule 
given in Eq. ([2]). Since the parameters a,; are infinitesi- 
mal, we have 7(01,02) = 2 and the pair of variables uji 
and fii transforms autonomously (i.e. not influenced by 
the variables ai), identically to the infection rate and re- 
covery rate of a symmetric contact process. As can be 
seen above, the asymmetry parameter o^ transforms in 
the same way as the length li of effective links given in 
Eq. ([To]). The only difference is in that the original pa- 
rameters Oj have not necessarily the same sign, while we 
have h = 1. The renormalized asymmetry parameter o; 
can thus be written as a sum of original variables Oi with 
li terms which are selected by the procedure indepen- 
dently of the parameters o; . 

The above renormalization scheme allows for two dif- 
ferent scenarios of the evolution of the asymmetry pa- 
rameters depending on the initial distribution of rates. 
There may be initial distributions biased to the left or to 
the right for which the asymmetry parameter ai averaged 
over sites tends to plus or minus infinity, respectively. Or 
it may tend to zero as the fixed-point of the transforma- 
tion is approached. In the latter case, the system is sta- 
tistically symmetric, although it is locally asymmetric. 
A sufficient condition for this is that the distribution of 
original infection rates is invariant under the interchange 
of A and k, i.e. 



P(A,k) = P(k, A). 



(21) 
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This property will then be preserved by the renormal- 
ization. It is, however, not a necessary condition for 
the system being statistically symmetric but owing to 
the approximative nature of the renormalization far from 
the asymptotical region and that even the homogeneous 
model is non-integrable it is not possible to find a general 
condition in terms of the distribution of the initial rates. 

On the basis of the above considerations, if the initial 
distribution of infection rates is biased so that the mean 
asymmetry is growing in magnitude, it must scale close 
to the fixed-point as the length: 



(22) 



But if the system is initially statistically symmetric, i.e. 
Q!(^) — > as the fixed-point is approached, then, accord- 
ing to the central limit theorem, the square of the vari- 
ance of a scales as the length: 



(23) 



Using now the results of the symmetric model quoted in 
the previous section, we obtain for the scaling of typical 
asymmetry parameter at the critical point A = Aq 



15(^)1 ^^''//2^ if 

5(0 ^ , otherwise. 



(24) 
(25) 



In the sub-critical Griffiths (A < Aq) phase, we obtain 

if 5->o (26) 

<5(C) ~ otherwise, (27) 
while in the super-critical Griffiths phase (A > Aq) 

|5(0I ~ (In^)'/', if (28) 
5(0 ^ In^, otherwise. (29) 



IV. THE STATISTICALLY SYMMETRIC 
MODEL 

Next, we consider finite local asymmetry parameters 
but require that the distribution of at is not too broad 
so that its moments are finite. As we have seen above, 
the asymmetry parameter scales differently depending on 
whether the system is statistically symmetric or not. Let 
us consider first, the former case, which is guaranteed if 
Eq. (EH) holds. 

Assume temporarily that the initial asymmetry is finite 
but weak (/ ^ 1) so that, at least within a finite length 
scale ^0 {I) I the system can be regarded as approximately 
symmetric. As the model is renormalized, the initially 
uncorrelated rates Aj and Ki become correlated since they 
are transformed simultaneously and on longer effective 
links both A and k are typically smaller than on shorter 
ones. If the system is critical, the typical asymmetry 
parameter starts to grow as \a\ = | InA — Inrej ^ I^'^f^^. 



But the typical value and the width of the distribution 
of logarithmic rates themselves are increasing faster, see 
Eq. ([TT|). since df/2 < 1/2. Consequently, if a combined 
infection rate w,; is much smaller than fi2 in a cluster 
elimination step, then typically both Ai and Ki are much 
smaller than fi2- In other words, a complete separation 
of niveous is realized. Furthermore, taking the logarithm 
of Eq. (|19p . we can sec that the term ln/(ai,a2), which 
is at most 0{S,'^f^^), is negligibly small compared to the 
other terms on the r.h.s., which are 0{£,^/^). Therefore 
the SDRG transformation can be continued beyond the 
scale S,oiI) and as the fixed-point is approached (^ — )■ 
oo), the relative asymmetry parameter tends to zero with 
probability one, 



In A,; 



and the system transforms asymptotically like the sym- 
metric one. But this is not restricted to the critical point. 
Following the above reasoning and using Eq. ()26|) . we 
obtain a vanishing relative asymmetry parameter in the 
sub-critical Griffiths phase. 



In Xi 



0, 



(A < Ao) (31) 



as well as, using Eq. ([25)) . in the super-critical Griffiths 
phase: 



\n{\i/ki) 
In A,: 



(lnC)-^/2->0. (A>Ao) 



(32) 



So, we conclude that, in the statistically symmetric 
model, the local asymmetry becomes irrelevant in the 
fixed-point of the SDRG transformation, and, conse- 
quently, the large scale properties are identical to that 
of the symmetric model. 

It is, however, difficult to answer whether for any type 
of initial distributions of transition rates the large scale 
properties of the system are described by the fixed-point 
of the SDRG transformation. For contributions to the 
unresolved question whether any weak initial disorder re- 
sults in logarithmic critical scaling described by the infi- 
nite randomness fixed-point in the symmetric model we 
refer the reader to Refs. [l^ [l^ [13, [2l| . In the asymmet- 
ric contact process this problem is made more difficult by 
introducing the extra parameter a^. The numerical inves- 
tigation of the SDRG scheme shows that not only weak 
but even relatively strongly asymmetric initial disorder, 
such as independent uniform distributions of the rates A^ , 
Ki and /i;, drives the system to the infinite randomness 
fixed-point. 



A. The site-symmetric model 

We make a digression here on a special model stud- 
ied by Monte Carlo simulations [l^, where the infec- 
tion spreads from each site in both directions with equal 
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rates, i.e. Xi = k^-i for all i. As opposed to the earlier 
assumption, the initial asymmetry parameters are now 
correlated variables and it is easy to see that the sum of 
Ui over subsequent sites is independent from the number 
of terms: 



y^Qi = ln(Ai/K„)- 



(33) 



i=l 



Thus, one expects here a slower increase of the asymme- 
try parameter under renormalization than for the uncor- 
related disorder. Assume for the sake of simplicity that 
the initial disorder is strong enough so that the asymp- 
totic forms of the renormalization rules can be used. The 
rcnormalizcd asymmetry parameter of an arbitrary effec- 
tive link is the sum of initial variables of links which 
have been incorporated into the effective link via cluster 
eliminations. These links, when regarded in the origi- 
nal system, form a disconnected set due to the eventual 
cluster merging steps occurring during the renormaliza- 
tion. Each connected part of length n of this set gives 
a contribution given in Eq. (|33p independently of n to 
d. So, the number of independent terms in a is given 
just by the number of connected components of the set 
of links rather than the total number of links. It is easy 
to obtain that the number of connected components q, 
which is initially 1 on each link, transforms in a cluster 
elimination step in the way: 



Cl + C2 — 1, if 712 = 1 
Ci + C2, if 7l2 > 1 



(34) 



This is not much different from the transformation of 
the variable li given in Eq. (|10p and, taking into account 
that the mass of clusters is increasing, the probability 
of finding a cluster with n = 1 goes to zero as ^ — >■ oo. 
So, the transformation rules of n and c are asymptotically 
identical and, consequently, the scaling of the asymmetry 
parameter follows the same laws as in the uncorrelated 
case given in Eqs. (|30ll32p . and at most the pref actors 
arc reduced. 



V. THE BIASED MODEL 



the forward rate A is typically small compared to /i, i.e. 
A//i — >■ 0. In this case, we have a complete separation of 
niveous in cluster elimination steps and the variables uj 
and fi transform autonomously as 



LJ1UJ2 



LJ2 



(35) 



in cluster elimination and merging steps, respectively, 
whereas the asymmetry parameter in the former as a ~ 
fli + 12- The rcnormalizcd parameters under these trans- 
formations will scale asymptotically as 



In k 



(36) 



and flow to a line of fixed points parameterized by a non- 
universal dynamical exponent z that depends on the dis- 
tribution of initial rates [l9j . This line can be interpreted 
as the subcritical Griffiths phase, where the dynamics are 
described by power laws with distribution-dependent ex- 
ponents. 

As in the Griffiths phase of the symmetric model, a fi- 
nite length scale ^1 above which the rcnormalizcd recov- 
ery rates are typically larger than the renormalized com- 
bined infection rates can be defined. This length scale can 
be interpreted as the characteristic size of locally super- 
critical clusters. Beyond this scale, ^ ^ ^1, practically no 
cluster merging step occurs in the SDRG procedure. In 
the biased model at this length scale ^ ^ i the backward 
rates k are comparable with fx but the forward rates A 
are still larger than p,. Similar to ^1, one can define a sec- 
ond length scale C2(> Ci) above which the recovery rates 
typically exceed even the forward rates. An interpreta- 
tion of ^2 can be given if the process starts from a single 
active site. In this case, the activity spreads to the right 
typically up to a finite distance ^2, where it is trapped 
until becoming extinct. Thus, the observables usually 
measured in simulations, such as the survival probability, 
are expected to have the same time dependence beyond 
a time scale corresponding to ^2 as in the sub-critical 
Griffiths phase of the symmetric model. 



In the remaining part of the paper wc will discuss the 
case when the spreading of activity is biased to either di- 
rection, for example, if > for all i, and will be inter- 
ested in the possible fixed-points of the SDRG transfor- 
mation. We have seen for infinitesimal asymmetry that 
q; — ^ 00 in the biased model. For finite asymmetry, we 
will therefore have typically Xi » ki on large scales and, 
consequently, cui ~ ki and f {0,1,0,2) ~ 1. 

A. Griffiths phase 

Assume now that the control parameter A is small 
enough, such that the system is sub-critical and even 



B. The lower phase transition and the 
weak-survival phase 

When the control parameter A is increased, both ^1 
and ^2 increase and there must be a point A = Ai, where 
^2 diverges but ^1 is still finite. This point corresponds 
to a transition to a phase where the activity spreads to 
the right without limits in a finite fraction of random, 
and the survival probability tends to a positive limit as 
t — > 00. At this point and above, p, will not exceed A 
asymptotically, therefore there is only a partial separa- 
tion of niveous. As aforementioned, we obtain then no 
correct information on how the activity spreads right- 
wards by the SDRG method and it is therefore unable to 
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predict the properties of this phase transition. The fact 
that the niveou corresponding to A does not separate 
suggests that the transition is, at least, not of activated 
type characterized by logarithmic time-dependence of ob- 
servables. Nevertheless, the lowest niveou separates well 
from the other ones even if A > Ai and in the asymp- 
totic transformation of the corresponding rate k and that 
of fj, the 'incorrect' rate A does not play a role. 

Although we cannot infer the properties of the phase 
transition, the spreading of the front of the activity 
started from a single active seed can be related to the dis- 
tribution of the renormalized recovery rates. If the sys- 
tem is renormalized well beyond the scale ^i, the renor- 
malized backward rates are negligibly small and it be- 
haves as a totally asymmetric contact process, where the 
recovery rates have a broad distribution with an algebraic 
tail: 

p<(Ai)~//^ ^i^o. (37) 

Assume now that the process was started from a single 
active seed at a;(i = 0) = and at time t there are still 
active sites. If the leftmost active site, the position of 
which is denoted by x{t), recovers, it will be practically 
not reinfected and x is shifted to the closest active site 
on its right. The expected value of the time of recovery is 
1/ fj-i, therefore the mean time ti.„+i needed for x to shift 
from site 1 to site n -I- 1 is at most the sum of recovery 
times: 

n 

<^Ai,r'- (38) 

i=l 

For large n, the asymptotical behavior of this sum is 
governed by the exponent z. z < 1, ^^^^ ^ n, 
whereas if z > 1, "^27=1 1^7^ This yields the follow- 

ing lower bound on the asymptotical displacement of the 
leftmost active site: 

x{t) > 0(t) if z <1 
x{t) > 0{t^/') if z>l. (39) 

Obviously, the position of the rightmost active site (the 
front of the activity) must move at least as fast as x{t) 
in samples which are surviving up to time t. If A > 
Ai, the forward infection rates are typically larger than 
the recovery rates in the renormalized model, therefore 
the density of active sites between x{t) and the front is 
expected to be finite, similar to the super-critical phase 
of the symmetric model. Consequently, if the leftmost 
active site recovers, the expected shift of x{t) is finite 
and the relation in pop holds as a proportionality. 

So, the active cluster moves at least as a power of 
time, which is different from the behavior of the unbi- 
ased model. If z > 1, which is realized by a weak initial 
bias, the above reasoning suggests the creeping motion 
of the critical active cluster with an asymptotically van- 
ishing velocity. This is in contrast to the homogeneous, 
biased model where the active cluster moves with a finite 
velocity. 



C. The upper phase transition and the 
strong-survival phase 

If the control parameter A is further increased in the 
weak-survival phase (A > Ai), the length scale ^i, as well 
as the dynamical exponent z increase and, at a second 
phase transition point A = A2, they diverge. Physically, 
this means that for A > A2, the activity, when started 
from a single active seed spreads also leftwards without 
limits. In the transition point, the variables w ~ k and fi 
are dual to each other and the fixed point of the SDRG 
transformation is an infinite randomness fixed-point with 
the asymptotic scaling of parameters: 

\\nk\^\\nfi\^e^\ (A = A2) (40) 

Whereas, above A2, we obtain a line of fixed points pa- 
rameterized by a distribution-dependent dynamical ex- 
ponent z', where the parameters scale as 

^^r"', |ln/2|^e (A>A2) (41) 

The phase transition at A = A2 can be also detected 
in the change of the finite-size scaling of the lowest non- 
trivial eigenvalue of Q in a finite, open system of size L, 
which is (asymptotically) correctly treated by the SDRG 
method. The inverse of this quantity that gives the time 
scale of reaching the absorbing state from the fully active 
one scales in the following way: 

^_^max(l,^)^ Ai<A<A2 

lnT-Li/2, A = A2 

In T - L. A > A2 (42) 

Note that in a finite but periodic system, the finite-size 
scaling of the lowest gap of Q is different from this. In 
that case, the extinction time behaves as In r ^ L in the 
entire range A > Ai, since the activity can 'go around' 
rightwards and it is irrelevant whether the spreading of 
infection leftwards is blocked or not. This is consistent 
with the SDRG treatment. When the periodic system 
is renormalized up to two effective sites then in the last 
decimation step the niveou which was the lowest (non- 
trivial) one till that point and which could indicate the 
phase transition at A2 is lost. 

As we mentioned above, in the biased model, the renor- 
malized forward rates will be much greater than the back- 
ward ones. So, it is plausible to assume that the behavior 
of the system in the range A > Ai is well described by a 
simply tractable idealized model where the forward rates 
are infinitely large Ai = 00. If this process is started 
from a single active site, the activity will immediately 
spread to the right without limits, and all sites on the 
right hand side of the leftmost active site will be active 
since once a site recovers, its left hand side neighbor will 
instantly reinfect it. It is easy to see that the position 
x{t) of the leftmost active site is a random walk with 
jump rate Pi = fii and qi ~ Hi-i to the right and to 
the left, respectively if x{t) = i. The basic properties of 
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the one-dimensional random walk in a random environ- 
ment are exactly known, for a review see e.g. Ref. p^ . 
Varying the control parameter of the problem, 



A 



RW 



\nqi - In Pi 



(43) 



the expected displacement x{t) in typical environments 
has different asymptotics: 



x{t) - t^/' 
\xit)\r^i\ntr 
-x{t) - i^/" 



Arw < 
Arw = 
Arw > 0, 



(44) 



1 



where z is the positive root of the equation {qi/pi)^/ 
for Arw < 0. 

So, in the range Ai < A < A2, the position of the 
leftmost active site is expected to creep rightwards, at 
A = A2 to perform a Sinai walk (23| with a zero aver- 
age displacement and fluctuations in the order of (Ini)^, 
finally, if A > A2, to creep leftwards. These results are 
consistent with the finite-size scaling of the extinction 
time given in Eq. (j42p and are in line with the behavior 
of the homogeneous, biased contact process, where x{t) 
is a random walk in a homogeneous environment and the 
extinction time scales at the second transition point as 



24 



The phase transition at A = A2 can also be detected in 
the behavior of the survival probability in a semi-infinite 
system, when the process is started from a single ac- 
tive site (site 0) and infection spreading from this site 
is not possible to the right, i.e. Aq = 0. This quantity 
can be calculated in the idealized model, where the site 
on the right hand side of the origin of the random walk 
will be an absorbing site. The probability P{x) averaged 
over the random environments (i.e. sets of transition 
rates) that the walker visits a site in a distance x from 
the origin before it hits the absorbing site is known to 
be 0(6-™"*^*-^), 0(a;-i/2) and constant for Arw < 0, 
Arw = and Arw > 0, respectively, see e.g. Ref. [25l |. 
Using these results, the asymptotical time-dependence of 
the survival probability can be derived. If A < A2, the 
activity is localized in the vicinity of the origin but, due 
to the algebraic tail of the distribution of effective recov- 
ery rates in Eq. p7p . the survival probability averaged 
over disorder decays algebraically, as well: 



Psurv(t)^i-'/^ (A<A2) 



(45) 



At the second critical point, A = A2, using the relation 
between time and length scale of the Sinai walk Int ~ 
which can be gathered also from the fixed-point 
solution of the SDRG transformation in Eq. (|40|. we 
obtain for the survival probability: 



-Psurv(<) - (InO' 



(A = A2) 



(46) 



This form is identical to the surface critical behavior of 
Psurv {t) hi the unbiased model . An important differ- 
ence is, however, that the set of active sites is compact 



in the biased model, i.e. has a finite density, while in the 
symmetric model it is a fractal of dimension df. 

Finally, for A > A2, the survival probability tends to 
a positive limit that depends on the initial distribution 
of rates when t ^ 00. 



VI. SUMMARY AND OUTLOOK 



We have applied an SDRG method to the asymmetric 
contact process where the infection rates to the left and 
to the right, as well as the recovery rates are indepen- 
dent random variables. We have shown that the local 
random asymmetry in the infection spreading is irrele- 
vant on large scales if the model is globally (statistically) 
symmetric, and the critical and off-critical behavior is 
identical to that of the disordered, locally symmetric con- 
tact process. If the model is globally biased, the SDRG 
transformation predicts two distinct phase transitions. 
The lower one is related to the transmission of infection 
in the direction of the bias. Unlike this transition point, 
which is out of the range of validity of the method, the 
upper critical point, which is related to the spreading of 
activity against the direction of the bias, is controlled by 
an infinite randomness fixed point. 

It would be desirable to check the predictions of the 
SDRG method by Monte Carlo simulations. An intrigu- 
ing question which is left open by the SDRG analysis is 
nature of the lower phase transition. This is predicted 
to be neither of activated type with logarithmic critical 
scaling nor to be in the directed percolation class at least 
for weak asymmetry when the critical cluster is creeping. 
The present work is restricted to d = 1, although the 
critical behavior of the symmetric contact process in the 
presence of disorder is known to governed by an infinite 
randomness fixed-point also in dimensions d > \ (26l428| . 
The question of stability of the fixed-point against asym- 
metry in higher dimensions remains a challenging open 
problem. 
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Appendix A: The spectrum of the rate matrix 
1. Two-site block 

The rate matrix of a two-site block shown in Fig. [T]is 
given as: 



where 



a = Ai + Ki + 2(/ii + ^2) 

b = AiKi + ^iyU2 + (/-ii + M2)(Aii + M2 + ^1 + ^i) 

C = AilyU2(Ail + + Ai + Ki). (A3) 



Q12 = 

















-Ml - ^1 





Ai 


tJ-2 





-M2 - Ki 







M2 


Ml 


-Ml - M2 



(AI) 



The three non-trivial eigenvalues are roots of the cubic 
equation: 



Three-site block 



The rate matrix of the three-site block shown in Fig. 



Q123 — 



+ be 


+ c = 0, 




(A2) 

1 


[2] has a 


block diagonal form: 


























M2 


— Ki — Ai — M2 





Kl 





A2 














-Ai 


Ai 




















Ml 


-A2- 


Ml 








A2 














-K2 


K2 




















M2 


-Kl - Ml 

























— Al — K2 


Al + K2 


Vo 

















M2 


-M2 / 



(A4) 



and has the following eigenvalues: 

ei = £2 = 0, 

|e3,4| = ^(M2 + Al + A2) X 



|e5,6| = 2(^2 + Kl + K2) X 



(A5) 



X (1 T V1-4AiA2/(m2 + Ai+A2)2) (A6) 



X (1 T -4kiK2/(m2 + Ki +^2)2) (A7) 
le?! = M2 + Al -f K2, lesi = M2 + A2 + Kl. (A8) 
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